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ABSTRACT 

Aims. Massive stars form in clusters, and they are often found in different evolutionary stages located close to each other. To understand 
evolutionary and environmental effects during the formation of high-mass stars, we observed three regions of massive star formation 
at different evolutionary stages that reside in the same natal molecular cloud. 

Methods. The three regions S255IR, S255N and S255S were observed at 1.3 mm with the Submillimeter Array (SMA) and follow- 
up short spacing information was obtained with the IRAM 30m telescope. Near infrared (NIR) H + A^-band spectra and continuum 
observations were taken for S255IR with VLT-SINFONI to study the different stellar populations in this region. 
Results. The combination of millimeter (mm) and near infrared data allow us to characterize different stellar populations within the 
young forming cluster in detail. While we find multiple mm continuum sources toward all regions, their outflow, disk and chemical 
properties vary considerably. The most evolved source S255IR exhibits a collimated bipolar outflow visible in CO and H2 emission, 
the outflows from the youngest region S255S are still small and rather confined in the regions of the mm continuum peaks. Also the 
chemistry toward S255IR is most evolved exhibiting strong emission from complex molecules, while much fewer molecular lines are 
detected in S255N, and in S255S we detect only CO isotopologues and SO lines. Also, rotational structures are found toward S255N 
and S255IR. Furthermore, a comparison of the NIR SINFONI and mm data from S255IR clearly reveal two different (proto) stellar 
populations with an estimated age difference of approximately 1 Myr. 

Conclusions. A multi-wavelength spectroscopy and mapping study reveals different evolutionary phases of the star formation regions. 
We propose the triggered outside-in collapse star formation scenario for the bigger picture and the fragmentation scenario for S255IR. 

Key words, stars: formation - ISM: jets and outflows - ISM: molecules - stars: early-type - Hertzsprung-Russell (HR) and C-M 
diagrams - stars: individual: S255IR, S255N, S255S 

^ ■ 1. Introduction The central region S255I R: The central IRA S source with 

H ' a luminosity of 5 x 10 4 L (Minier et al. 2005) harbors three 
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UCHll regions region (Snell& Bally 1986) that is associated 

Massive stars (M > 8M ) are one of the paramount com- with Class II CH3OH and H 2 maser emission (IGoddi et alj 

ponents in the evolution of the universe, yet their formation 2007), which indicate the presence of massive young stellar 

is significantly less well understood, than that of their low- objects. The luminosity would be 2 x 10 4 L if we apply the 

mass counterparts. S255IR is a famous m assive star forma - new distance of 1.59 kpc. The region hosts a cluster of low- 

tion complex at a distance of 1.59^ ^ kpc (iRygl et al.ll2010h . mass sources that is surrounded by a shocked bubble of H 2 gas 

embraced by the Sharpless regions S255 and S257 that are (lOiha et al.ll2006t iMiralles et al.lfl997h . iTamura et all (Il99ll) re- 

alrea dy evolved Hll regions. The SCUBA 850 fim observa- solved the central near infrared sou rce into two compac t sources 



alrea dy evolved Mil regions. 1 ne SLUM S3U /im observa- soivea me central near inrrarea source into two compact sources, 
tion dDi Francesco et all 120081: iKlein et alJ | 2005l) shows three MRS 1 and NIRS 3. Furthermore. iMinier et all d2007l) reported 
main continuum sources: G192.60-MM2 (Mi nier et al.l 120051) HCO + infall signatures toward this region. In the mid-infrared, 



signatures 

lies in the center region S25 5IR, and two additi onal mm con- lLongmore etafl d2006l) resolved a massive proto-binary which 
tinuu m peaks G192 .60-MM1 (IMinier et al.ll2005l) and G192.60- coincide with NIRS 1 and NIRS 3 (Figure©. And NIRS 1 has 
MM3 ( IMinier et al.l2005l) toward the northern region S255N and been identified as a massive disk candidate by near-infrared po- 
southern region S255S, respectively (Figure[TJ. 
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lariza tion observations (Jiang et al. l2008l see also lSim pson et al. 
(2009) for a different interpretation). 

The northern region S255N: This region with 10 5 L 
dMinier et al .1120051) also hosts a UCHll region G192.584-0.041 
dKurtz et al.lll994h as well as Class I CH^OH and H?Q maser 
emission (Kurtz et al. 20 041 ICveanowski et al.l l2007). which in- 
dicate the presence of the massive young stellar objects. The 
luminosity would be 4 x 10 4 L if we apply the new distance 
of 1.59 kpc. The recent Submillimeter Array observations to- 
ward that region at a spatial resolution of 3.6" resolved three 
mm sources and strong molecular line emissi on, and none of the 
sourc es is associated with NIR point source (Cvganowsk i et alj 
120071) . Outflow activity is evidenced by various tracers from 
shocked H? emission via cm continuum a nd SiO emission 
dMiralles et alj|1997t ICyganowski et al]|2007l) . and global infall 
was reported bv lMinier et al.l d2007l) . 

The southern region S255S: This sub-source is the least stud- 
ied so far. It exhibits strong mm continuum emission (Figure [TJ 
but no other sign of active star formation yet such as near in- 
frared and mid infrared emissions. Therefore, it is p roposed to be 
in a p re-stellar phase of the evolutionary sequence dMinier et alJ 
120071) . 

While G192.60-MM2 in S255IR still shows signs of active 
star formation, it is associated with a NIR cluster and appears to 
be the most evolved one of the three region. G192.60-MM1 in 
S255N is of similar luminosity, and G192.60-MM3 is S255S is a 
High-Mass Starless Core candidate. Therefore, S255 complex is 
the ideal candidate source to simultaneously investigate several 
sites of massive star formation at different evolutionary stages 
within the same larger-scale environment. 

2. Observations and Data reductions 

2. 1 . Submillimeter Array observations 

The S255 complex was observed with three fields with the 
Submillimeter Array (SMA) on November 3rd 2008 in the 
compact configuration with seven antennas and on February 
8th, February 17th 2009 in the extended configuration with 
eight antennas, and February 13th in the extended con- 
figuration with seven antennas. The phase centers of the 
fields, which are known as S255IR, S255N and S255S, 
were R.A. 06hl2m54.019s Dec. +17°59'23.10" (J2000.0), 
RA. 06hl2m53.669s Dec. +18°00'26.90" (J2000.0) and R.A. 
06hl2m56.58s Dec. +17°58'32.80" (J2000.0), respectively. The 
SMA has two spectral sidebands, both 2 GHz wide and separated 
by 10 GHz. The receivers were tuned to 230.538 GHz in the up- 
per sideband (v/ sr =10 km s ) with a maximum spectral resolu- 
tion of 0.53 km s . The weather of February 8th and February 
13th was mediocre with zenith opacities t(225 GHz) larger than 
0.2 measured by the Caltech Submillimeter Observatory (CSO). 
But the weather of November 3rd 2008 and February 17th 2009 
was good with zenith opacities r (225 GHz) around 0.1. So we 
only used the data observed on these two days. For the com- 
pact configuration on November 3rd 2008, bandpass was de- 
rived from the quasar 3c454.3 observations. Phase and ampli- 
tude were calibrated with regularly interleaved observations of 
the quasar 0530+135 (11.4° from the source). The flux calibra- 
tion was derived from Uranus observations, and the flux scale is 
estimated to be accurate within 20%. For the extended config- 
uration on February 17th 2009, bandpass was derived from the 
quasar 3c273 observations. Phase and amplitude were calibrated 
with regularly interleaved observations of the quasar 0530+135. 
Because of the lack of flux calibrator observations, the flux was 



estimated by the SMA calibrator database for the gain calibra- 
tor, to be accurate within 20%. We merged the two configura- 
tion data set together, applied different robust parameters for the 
continuum and line data, and got the synthesized beam sizes be- 
tween 1.4" xl.l" (PA 86°) and 1.9" x 1.6" (PA 77°), respec- 
tively. The 3<x rms of 1 .3 mm continuum image is ~ 4.5 mJy and 
the 3cr rms of the line data is 0.14 Jy/beam at 2 km s _1 spec- 
tral resolution. The flagging and calibr ation was done with the 
IDL superset MIR ( Sc oville et al.l [19931) which was originally 
developed for the Owens Vally Radio Observatory and adapted 
for the SMAB The imaging and data analysis were conducted in 
MIRIAD (ISaultet al.ll 1991 . 

2.2. Short spacing from the IRAM 30 m 

To complement the CO(2 - 1 ) observations with the missing short 
spacing information and to investigate the large-scale general 
outflow properties to find the connection of the three regions, we 
observed them with the HERA array at the IRAM 30 m telescope 
on November 10th 2009. The 12 CO(2 - 1) line at 230.5 GHz and 
I3 CO(2 - 1) line at 220.4 GHz were observed in the on-the-fly 
mode. We mapped the whole region with a size of 5'x 3' cen- 
tered at R.A. 06hl2m54.02sDec. +17°59'23.10" (J2000.0). The 
sampling interval was 3.5", a bit better than Nyquist-sampling to 
minimize beam-smearing effects. The region was scanned two 
times in the declination and right ascension direction, respec- 
tively, in order to reduce effects caused by the scanning process. 
The spectra were calibrated with CLASS which is part of the 
GILDAS software package, then the declination and right as- 
cension scans were combined with the plait algorithm in GREG 
which is another component of the GILDAS package. The I2 CO 
data has a beam size of 11", and the rms noise level of the cor- 
rected T m b is around 0.23 K at 0.6 km s _1 spectral resolution, the 
I3 CO has the same beam size but the rms noise is 0.12 K at 0.6 
km s _1 spectral resolution. 

After reducing the 30m data separately, single dish I2 CO and 
I3 CO data were converted to visibilities and then combined with 
the SMA data using MIRIAD package task UVMODEL. With 
different uv-range selections, the synthesized beam of the com- 
bined data varies from 2.1" x 1.8" (PA -84°) to 11.1" x 10.7" 
(PA -74°). 

2.3. VLT-SINFONI integral field spectroscopy observations 

Near Infrared H- and A"-band observations we re performed us- 
ing the Integral Field in strument SINFONI (Eisenhaue r et al.l 
2003: lBonnet et al.l20 04) on UT4 (Yepun) of the VLT at Paranal, 
Chile. The observations of S255IR were performed in service 
mode on February 8th, February 12th, February 13th and March 
29th 2007. The non-AO mode of SINFONI, used in combina- 
tion with the setting providing the widest field of view (8" x 8") 
with a spatial resolution of 0.25" per slitlet. The typical seeing 
during the observations was 0.7" in A"-band. The H+K grating 
was used, covering these bands with a resolution of R=1500 in 
a single exposure. 

S255IR was observed with a detector integration time (DIT) 
of 30 seconds per pointing. The observations were centered on 
the central near infrared source NIRS 1 at coordinates: R.A. 
06hl2m53.85sDec. +17°59'23.71" (J2000.0). We observed this 
area with SINFONI using a raster pattern, covering every loca- 
tion in the cluster at least twice, resulting in an effective integra- 

1 The MIR cookbook by Chunhua Qi can be found at 
http://cfa-www.harvard.edu/~cqi/mircook.html 
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Table 2. SCUBA sub-millimeter continuum data. 



Source 


S peak 




Mass 


N„ 2 




Jy/beam 


Jy 


[M Q ] 


cirr 2 


S255IR 


7.81 


31 


372 


6.5 x 10 23 


S255N 


7.75 


30 


357 


6.4 x 10 23 


S255S 


2.33 


9 


110 


1.9 x 10 23 



tion time of 60 seconds per location in the field. The offset in the 
east-west direction was 4" (i.e. half the FOV of the detector) and 
the offset in the north-south direction was 6.75", results a total 
observation field of 70" x 70". A sky frame was taken every 3 
minutes using the same DIT as the science observations. The sky 
positions were chosen based on existing NIR imaging in order 
to avoid contamination. Immediately after every science obser- 
vation, a telluric standard star was observed, matching as close 
as possible the airmass of the object. 

The data were reduced using the SPRED (version 
1.37) software developed by the MPE S INFONI consortium 
(ISchreiber et alJ 12004 lAbuter et all l2006h . The procedure de- 
scribed by iDavies i d2007l) was applied to remove the OH line 
residuals. To calibrate the flux and remove the telluric absorp- 
tion lines, the extracted spectrum of one standard star in each OB 
(OBI: Hip064656, OB2: Hip048128, OB3: Hip032108, OB4: 
Hip0318 99) is used. The flu x calibration of the spectra uses the 
2MASS dCohen et al.ll2003h magnitude of the standard stars (see 
iBik etalJd20irI for the details of the data reduction). 



3. Results 

3. 1 . Millimeter continuum emission 

We averaged the apparently line-free parts of upper and lower 
sideband spectra of each region (presented in Figure |2]i and got 
the continuum images of the three regions shown in Figure [3] 

Assuming optical thin dust emission, we estimated the gas 
mass and column densit y of the con t inuum sour ces following 
the equ ations outlined in Hildebrand (1983) and Beuther et al.l 
(2005a). We assumed a dust temperature of 40 K, grain size of 
0.1 /mi, grain mass density of 3 g cm 3 , gas-to-dust ratio of 100 
and a grain emissi vity index of 2 (correspondi ng to k « 0.3 
for comparison of lOssenkopf & Henningl {[994)). The results 
are shown in Table Q] With the same set-up, we also calculated 
the gas mass and the column density of th e SCUBA 850 /mi 
dDi Francesco et aD 120081: iKlein et all 120051) continuum peaks 
shown in Table [2] To get the impression of how much flux we 
lost in the interferometer observations, we produced the SMA 
continuum map of the each regio n with the same beam size as the 
SCUBA 850 /mi map (14"xl4". |p"i Francesco etai1(l2008l» . We 
measured the continuum flux, calculated the gas mass and com- 
pare the results with the ones in Table |2j the mass we obtained 
from the SMA observations for S255IR, S255N and S255S is 
only 7.8%, 8.7% and 2% of the SCUBA 850 /im measurements, 
respectively. Since our interferometer observations are not sen- 
sitive to spatial scales > 24" (the shortest base line =8.5 kA), 
we filter out the smoothly distributed large-scale halo and left 
only the compact cores. Also, the fact that we filtered out more 
flux for the youngest region S255S indicates that at the earli- 
est evolutionary stages the gas is more smoothly distributed than 
during later stages where the collapse produces more centrally 
condensed structures. Nevertheless, all derived column densities 
are of the order or above t he proposed threshold for high-mass 
star formation of 1 g cm ~ 2 dKrumholz & M cKee 2008). 




6 h 12 m 55" 50" 
RA CJZCOO) 



6 h 12 m 55 B 50" 
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Fig. 1. SPITZER IRAC 4.5 /mi image (jChavama et al. 2008) 
overl a id with SCUBA 850 //m contours dDi Francesco et all 
2008: IKlein et al.ll2005h . In the upper panel, the grey scale is 
the SPITZER IRAC 4.5 //m image, the black contours are the 
NVSS 1.4 GHz emission and the green contours are the SCUBA 
850 /mi continuum. In the bottom left panel, the grey scale is 
the SPITZER MIPS 24 /mi image, the contours are the SCUBA 
850 /mi continuum, the three dashed circles mark the primary 
beam of our SMA observations. In the bottom right panel, the 
grey scale is the SPITZER MIPS 70 /mi image, the contours and 
the circles are the same as the ones in the left panel. The con- 
tour levels of the NVSS start at 10cr (lcr = 0.6 mJy beam -1 ) 
with a step of 5<x. For the SCUBA 850 /mi, the contour lev- 
els start at 10<x (lcr = 0.1 Jy beam -1 ) with a step of 10cr. The 
SPITZER/IRAC post-bed data processed with pipeline version 
SI 8.7.0 and the MIPS post-bed data processed with pipeline ver- 
sion S 17.2.0 have been downloaded from the SPITZER archive 
to create these images. 



S255IR The continuum image of S255IR is presented in the 
left panel of Figure [3] We resolved two main continuum peaks 
(i.e. S255IR-SMA1 and S255IR-SMA2) and one unresolved 
sub peak (i.e. S255IR-SMA3) in this region. The stronger peak 
named as S255IR-SMA 1 in the southeast co incides with the near 
infrared source NIRS 3 dTamura et alJ 19911) S255 IR-SMA1 co- 
incides with a UCHll region ( iSnell & Ballylll986l) generated by 
NIRS 3, and the inset in the S255IR panel in Figure [3] shows 
that it is also associated with Class II CHjOH and H2O maser 
emissions dGoddi et alj|2007t iMinier et alj feoOO, 2003- An un- 
resolved peak S255I R-SMA3 that coinci des with the near in- 
frared source NIRS 1 dTamura et alJI 199 ll) . which lies 2.4" west 
to S255IR-SMA1 and has been identifi ed as a massive d isk can- 
didate by NIR polarization observation (Jian g et al.l2008l) . These 
two continuum sources are both respectively coincident with the 
mid-infrared massive p roto-binary source 1 and 2 identified by 
lLongmore et al.l d2006l) . S255IR-SMA2, which has never been 
detected before, is not associated with any infrared source and 
likely to be in a very young source, it also shows only a few 
lines, which will be discussed in detail in Sec. 
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Table 1. Millimeter continuum peaks properties. 



S m 1 rpp 




R.A. 


Dec. 


^ peak 




Mass 






N„ 2 




(J2000.0) 


(12000 s ) 


m Tv/hpam 

1 i i.i y j i-.>\- ill 1 1 


mjy 


r/w„i 

l Iyl oi 




cm - 


J>Z351K-5>MA1 


06 


12:54.01 


+ 1 l.jy.li.l 


93 


1 "7 T 
1/1 


12 


3.3 


X 


10 24 


S255IR-SMA2 


06 


12:53.77 


+ 17:59:26.1 


52 


150 


10 


3.0 


X 


10 24 


S255IR-SMA3 


06 


12:53.88 


+ 17:59:23.7 


30 


30 


2 


1.7 
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10 24 


S255N-SMA1 


06 


12:53.71 


+18:00:27.3 


91 


327 


22 


5.5 
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10 24 


S255N-SMA2 
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Fig. 3. The SMA 1.3 mm continuum map of S255IR, S255N and S255S. contour levels start at 5cr with a step of 5<x in S255IR 
(lcr = 1.7 mJy/beam) and S255N (lcr = 1.6 mJy/beam) image and 3cr in S255S (lcr = 0.9 mJy/beam) image. The dotted contours 
are the negative features due to the missing flux. The synthesis beam is sh own at the bottom le ft of each figure. In the S255IR panel, 
the asterisk marks the position of NIRS 3 and the right one marks NIRS 1 dTamura et alj 199lh . The inset is a zoom in the continuum 
image of the inner region of S255IR-SMA1 and the contour levels start at 5cr with a step of 10<r. The triangles are the water maser s 
(Go ddi etal] l2007) and the two open diamond mark the positions of two 6.7 GHz methanol masers detected by IXu et alJ (12009). 
In S255N panel, the crosses mark the position of the Class I 44 GHz (7o - 6i) A + methanol masers detected by iKurtz et al.l ([2004) 
and the triangle is the water maser dCyganowski et al]|2007t) . The (0,0) point in each panel from left to right is R.A. 06hl2m54.019s 
Dec. 17°59'23.10" (J2000.0), R.A. 06hl2m53.669s Dec. +18°00'26.90" (J2000.0) and R.A. 06hl2m56.58s Dec. +17°58'32.80" 
(J2000.0), respectively. 



S255N The middle panel of Figure [3] is the continuum image 
of S255N. The strongest peak named as S255N-S MA1 which 
coinc ides with the UCHll region G192.584-0.041 dKurtz et alJ 
119941) is associate with Class I 44 GHz (7p - 6i) A + methanol 
masers which are detected by IKurtz et all d2004l) and water 
maser (Cvganowski et al 112007b . All the methanol masers which 
are marked with crosses in the S255N panel in Figure[3]are dis- 
tributed along the direction of the outflow (Figure [TO] see Sec. 
13. 3I >. S255N-SMA1 is also assoc iated with centimeter contin- 
uum emission (ICvganowski et alj 2007) which shows an elonga- 
tion aligned with the outflow (FigureQJJl left panel). However, all 
three continuum sources in S255N have no corresponding near 
infrared point sources. 

S255S The right panel of Figure [3] presents the continuum im- 
age of S255S. There are no near infrared sources or mid infrared 
peaks (SPITZER MIPS 24 fim and 70 //m, Figure [Qi that are 
associated with these two sources. Meanwhile, the gas mass is 
only 2% of the SCUBA 850 yum measurements, which indicates 
that the gas in this region is at an early stage and smoothly dis- 



tributed. All these features indicate that S255S is an extremely 
young region. 

3.2. Spectral line emission 

The observed spectra of lower sideband (LSB) and upper side- 
band (USB) of three region are shown in Figure|2] 

S255IR In S255IR, we detect 25 lines from 10 species. Besides 
three normal CO isotopologues, we also detect some sulfur- 
bearing species (SO, 34 S02, OCS, 13 CS) and some dense gas 
molecules which are usually used to trace hi gh mass hot cores, 
such as CH^OH, CH^CN and HCOOCH^ dNomura & Millarl 
12004 iBeuther et al]|2009h ISutton et al.lll985l) . All the lines we 
detect are with lower energy levels Ei ower /k between 5.3 to 568 
K (TableEJ. 

Figures |4] and present the integrated line images of all 
species (except 34 S02 which is blended with CH3OH and too 
weak for imaging, and 12 CO and I3 CO will be discussed in de- 
tails in I3.31 > with different transitions. Most lines show compact 
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Table 3. Observed lines in S255IR. 



2 
1 .5 



0.5 



Upper Sideband 
o CH 3 0H 



ocs cs : 



229.5 230 230.5 231 

Frequency [GHz] 

Fig. 2. Lower and Upper sideband spectral vector-averaged over 
all baselines with a resolution of 2 km s _1 per channel. In the 
Upper sideband image, the I2 CO line in S255IR and S255N pan- 
els was not fully plotted, which goes to 1 1 .6 Jy and 5.9 Jy respec- 
tively. 
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CH 3 OH(3 2 . 2 -4,, 4 )E 


28 


230.538 


12 CO(2-l) 


5.5 


231.061 


OCS(19-18) 


100 


231.221 


13 CS(5 -4 ) 


22 


231.281 


CH 3 OH(10 2 ,,-9 3 , 6 )A- 


98.8 


. Observed lines in S255N. 


V 


line Slower//! 


[GHz] 




[K] 


LSB 


219.560 


C ls O(2-l) 


5.3 


219.949 


SO(6 5 -5 4 ) 


24 


220.079 


CH 3 OH(8o, 8 -7 I , 6 )E 


85 


220.339 


13 CO(2-l) 


5.3 


220.743 


CHsCNC^-lU) 


65 


220.747 


CH 3 CN(12 -11 ) 


58 


USB 


229.759 


CH 3 OH(8_i,g-7o, 7 )E 


77 


230.538 


12 CO(2-l) 


5.5 


231.061 


OCS(19-18) 


100 


231.221 


13 CS(5o-4 ) 


22 



emission peaked at S255IR-SMA1, and only a few show emis- 
sion at S255IR-SMA2. Several CH 3 OH lines (CH 3 OH(8 g - 
7i, 6 )E, CH 3 OH(8-i, 8 - 7 ,v)E and CH 3 OH(3 2 , 2 - 4 1A )E, in 
Figure |4|i exhibit extended emission towards S255IR-SMA2. 
The CH30H(8_i 8 - 7o,7)E line even extends to the northwest 
of S255IR-SMA1 following the direction of the blue-shifted 
outflow of S255IR-SMA1 (Figure E see Sec. |3Jj. We sug- 
gest that particular the (8_i,8 - 7o ; 7)E line is due to the shock 
heating excited by the outflow, which could also be mixed up 
with some Class I methanol maser emissi on dSutton et all 2004: 
ISobolev et alj2005tlKdenskii et al.l2002tlSlvsh et al.l2002l) . SO 
line emission associated with S255IR-SMA1 is elongated along 
the outflow and could be affected by the outflow. 13 CS line emis- 
sion forms a shell around S255IR-SMA1 (Figure[5| simil ar shei l- 
like C 34 S emission has also been found by Beuthe r et al.l ([2009)). 



Only a few lines are associated with S255IR-SMA2, SO, C lg O 
(Figure|5]l and several methanol lines (Figure|4|i. S255IR-SMA2 
appears to be chemically younger than S2555IRmml. We ex- 
tracted C ls O spectra from these two continuum peaks. The spec- 
trum toward S255IR-SMA1 shows two peaks and can not be 
fitted by a single Gaussian profile, therefore, we calculated the 
FWZI (full width at zero intensity) of these two C ls O spec- 
tra, which is 12 km s" 1 and 6.8 km s _1 for S255IR-SMA1 and 
S255IR-SMA2, respectively. The difference of the line width 
also suggests that S255IR-SMA2 is at a younger evolutionary 
stage. 

S255N In S255N, we detected 10 lines from 6 species ( 12 CO, 
SO, CH3OH, CH 3 CN, OCS, 13 CS) and 2 additional CO isotopo- 
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Fig. 4. S255IR CH3OH line integrated intensity images with the SMA 1.3 mm continuum emission in the background. Contour 
levels start at 3<x with 2cr/level. For the top panels, the cr of the contours in each panel, from left to right, is 12, 15, 1 1 mJy beam -1 , 
respectively, and for the bottom panels is 12, 14, 15 mJy beam -1 , respectively. The dotted contours are the negative feature. The 
integral velocity ranges are shown at the bottom-left of each panel in km s . The synthesized beams are shown at the bottom-left 
of each panel. The (0,0) point in each panel is R.A. 06hl2m54.019s Dec. +17°59'23.10" (J2000.0). 



Table 5. Observed lines in S255S. 



V 


line 




[GHz] 




[K] 


LSB 


219.560 


C 1S 0(2-1) 


5.3 


219.949 


SO(6 5 -5 4 ) 


24 


220.339 


13 CO(2-l) 


5.3 


USB 


230.538 


12 CO(2-l) 


5.5 



logues ( 13 CO and C 18 0) with lower energy levels Ei OHVr /k be- 
tween 5.3 to 65 K (TableHJ. Figure [6] shows the integrated line 
images of all species (except three CO isotopologues which will 
be discussed in Sec. I3.3l and l3.4l >. All the lines show a peak as- 
sociated with S255N-SMA1, CH3OH and SO lines also show 
extended emission along the direction of the outflows (Figure 
[TO] see Sec. O- 

S255S In S255S, we detected 4 lines from two species ( 12 CO 
and SO) and 2 additional CO isotopologues ( 13 CO and C ls O) 



with lower energy levels Ei ower /k between 5.3 to 24 K (Table|5]l, 
which again indicates that S255S is an extremely young and cold 
region. We measure the single dish 13 CO line width toward the 
S255S SCUBA 850 //m continuum peak, and calculate the virial 
mass in this region, which is ~ 80 M Q . The virial mass is smaller 
than the mass we got from the SCUBA 850 fim measurement 
which implies that S255S region will likely collapse. 

From the single dish 13 CO observations, we got the v/ jr of 
the three regions, which are shown in Table [6] However, the 
v/ jr of the high-mass mm cores in S255IR and S255N based 
on the interferometer observations is different from the one we 
got from the single dish data. We extracted dense gas spectra, 
i.e. CH 3 OH for S255N-SMA1, CH 3 CN for S255IR-SMA1 and 
C ls O for S255IR-SMA2, and then a Gaussian profile was fitted 
to get the v/ sr of the cores. All the velocities are shown in Table 
[6] For S255S we can not get a proper v/ jr from the SMA obser- 
vations due to the missing flux problem, so only the v; sr we got 
from 30m observation is listed. 

3.3. The molecular outflows 

After combining the single dish 30m data with the interferome- 
ter data, we integrated the line-wing of 12 CO and produced the 
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Fig. 5. S255IR molecular line integrated intensity contour with the SMA 1.3 mm continuum emission in the background. All the 
contour levels start at 3 a. The contour step is 3 cr in SO (lcr = 22 mJy beam 1 ), HCOOH (lcr = 10 mJy beam 1 ), HCOOCH3 
(lcr= 12 mJy beanr 1 ), C ls O(lcr = 16 mJy beanr 1 ), CH 3 CN (lcr = 1 1 mJy beanr 1 ) and HNCO (lcr = 13 mJy beanr 1 ) images, 
and 1 cr in OCS (lcr = 13 mJy beam _1 )and 13 CS(lcr =11 mJy beam -1 ) image. The dotted contours are the negative features 
with the same contours as the positive ones ini each panel. The integral velocity ranges are shown at the bottom-left of each panel 
in km s _1 . The synthesized beams are shown at the bottom left of each panel. We only plotted the k = 0&1 line for CH3CN and 
(lOo.io - 9 0)9 ) line for HNCO. The (0,0) point in each panel is R.A. 06hl2m54.019s Dec. +17°59'23.10" (J2000.0). 



outflow image of each region. Figure Q shows the single dish 
only outflow map. While the northeast-southwest outflow which 
is associated with S255IR is clearly shown in Figure [7] the high 
velocity outflow is also found in S255N region, however, the di- 
rection of the outflow could not be well defined in single dish 
outflow map. In the youngest region S255S, the single dish out- 



flow map does not show much emission. We combine the SMA 
data and 30m data together to study the outflow properties. 



S255IR Figure [8] shows the combined SMA and 30m outflow 
images of S255IR. The velocity-integration regime of the blue- 
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Fig. 6. S255N molecular line integrated intensity contour with the SMA 1.3 mm continuum emission in the background. All the 
contour levels start at 3 cr. The contour step is lcr in CH30H(8o,s - 7i 6)A- (lcr = 9 mJy beam -1 ) and CH3CN(fc = 0&1, lcr = 14 
mJy beam -1 ) images, 2cr in 13 CS (lcr = 25 mJy beam -1 ) image, 3 cr in SO (lcr = 13 mJy beam -1 ) and CH30H(8i8 -7o/7)E images, 
and 0.5 cr in OCS (lcr = 10 mJy beam -1 ) image. The dotted contours are the negative features. The integral velocity ranges are 
shown at the bottom-left of each panel in km s -1 . The synthesized beams are shown at the bottom left of each panel. The (0,0) point 
in each panel is R.A. 06hl2m53.669s Dec. +18°00'26.90" (J2000.0). 



Table 6. The v; sr measured with different telescopes. 



Sources 


V;„- got from 30m 


V; sr got from SMA 




km s -1 


km s -1 


S255N-SMA1 


8.6 


10.4 


S255IR-SMA1 


7.7 


5.2 


S255IR-SMA2 


7.7 


9.3 


S255S 


6.6 





shifted part is [-40,0] km s -1 and [16,56] km s -1 for red-shifted 
part. We selected different uv-ranges resulting in different reso- 
lutions. Although at different resolution the structure of the out- 
flow changes a lot, the northeast-southw est (NE-SW) outflow is 
confirmed in all images (also reported bv lMiralles et al.l d!997l) ). 
In the bottom-left panel, the H2 jet-like sources (region (a) and 
(b) in Figure[l9] see Sec. [3T67TT > also follow the NE-SW outflow 
direction. The water masers which are marked with triangles in 
the inset of S255IR panel in Figure [3] also follow the direction 
of the outflow, which confirms that S255IR-SMA1 is the driving 
source of this NE-SW outflow. In all panels except the bottom- 
right two, the red-shifted part of the outflow bends a little toward 



northwest, which might be a signature of the precessing jets. In 
the S255IR panel in Figure [3] the water masers also follow the 
NE-SW direction. We believe that NIRS 3 is the driving source 
of this outflow. As we are applying different uv-range selection, 
two red-shifted outflow components reveal themselves out to 
the north and south of the continuum sources, respectively, and 
these two components are both only shown in the lower velocity 
regime which is offset from the v/ s , 10 to 20 km s -1 . Finally, in 
the 30m-only panel and Figure [7] the red-shifted outflow com- 
ponent to the north of the continuum sources becomes connected 
with the component in S255N. 

Figure [9] presents the position-velocity (pv) diagram of the 
molecular gas along the outflow. The diagram shows that the 
red- and blue-shifted outflow resembles the Hubble-law with in- 
creasing velocity at longer distance from the outflow center, and 
the blue-shifted ou tflow also rese mbles the jet-bow-shock gas 
entrainment model dArce et alj|200l . 

S255N Figure [TO] shows the combined SMA and 30m outflow 
images of S255N. The velocity-integration regime of the blue- 
shifted part is [-44,0] km s -1 and [16,48] km s -1 for the red- 
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Fig. 7. The single dish outflow map of the S255 complex. The 
blue-shifted outflow is shown in dashed contours with a velocity- 
integration regime of [-40,0] km s _1 , and the red one is shown 
in full contours with a velocity-integration regime of [16,56] km 
s _1 . The contours start at Act with a step of 5 cr{\cr-\.2 K km 
s _1 ). The three crosses mark the position of the three SCUBA 
850 jum continuum sources. 



shifted part. The northeast-southwest (NE-S W) outflow is shown 
clearly in the SMA and SMA combined with 30m panels but 
is hardly seen in the 30m-only panel. In the middle panel, the 
crosses mark the position of the Class I 44 GH z (7o - 61) A + 
methanol masers detected by Kurtz et al. (2004) and the triangle 
is the water maser (Cv ganowski et al]|2007l) . In the left panel, 
the back ground grey sca le is VLA Q-band 3.6 cm continuum 
dCvganowski et al] l2007) which shows a jet-like structure. All 
these features follow the NE-SW direction. Based on these fea- 
tures we confirm the direction of the main outflow as NE-SW. 
The red-shifted outflow to the north of S255N-SMA1 on the line 
(b) top left panel in Figure [10] might be part of the outflow cav- 
ity, but we can not exclude the possibility of multiple outflows. 
The red-shifted gas at the bottom part of line (b) seems to be 
associated with S255N-SMA3 (the middle panel of Figure [TOb, 
however, another blue-shifted feature at the southern part of the 
map shows up as the resolution changes. These two components 
both can also only be detected at the relatively lower velocity 
regime ~10 to 20 km s offset from the v/ s ,.. 

FigureQT]presents the position-velocity (pv) diagrams of the 
molecular gas along the outflows. In the top panel, the pv-plot 
cut follows the direction of the VLA jet-like emission (Fig.fTUl 
left panel, line a). The diagram shows that the high-velocity gas 
on the blue-shifted side remains very close to the outflow cen- 
ter and the red-shifted part does not show the Hubble-law sig- 
nature. But in the bottom panel, in which the pv-plot cut fol- 



lows the direction of the two elongate red-shifted emission (Fig. 
[TUl left panel, line b), the red-shifted northern side of the out- 
flow resembles the Hubble-law. However, the red-shifted south- 
ern part of the outflow seems to have nothing to do with our 
S255N-SMA1 nor S255N-SMA3. The different blue/red pv- 
characteristics may be explained by the blue 12 CO emission trac- 
ing predominantly the jet like component also visible in centime- 
ter continuum emission, whereas the red 12 CO emission traces 
mainly the cavity-like walls of the outflow or another outflow. 



S255S Figure Q~2] presents the combined SMA and 30m out- 
flow images of S255S. The velocity-integration regime of the 
blue-shifted part is [-9,1] km s _1 and [14,32] km s _1 for the 
red-shifted part. All the panels except the 30m-only one, show 
two blue-shifted components and one red-shifted component. 
Although the line wing emission is less pronounced than that 
for the other two regions, we clearly identify blue- and red- 
shifted 12 CO emission associated with both mm peaks. The blue- 
shifted part in the south of the map is only detected in the 1 km 
s~'>v> -6 km s _1 , and the red-shifted part in the north of the 
map is only detected in the v<16 km s _1 . Since the line wing 
emission does not extend to very high velocities, the outflow 
should not be oriented directly along the line of sight. Therefore, 
the small spatial extent of the outflow indicates as well the youth 
of this region. 
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Fig. 9. Position-velocity diagrams of S255IR for the 12 CO(2 - 1) 
SMA-only outflow observations with a velocity resolution of 1 
km s" 1 . The pv-diagram is in northeast-southwest direction with 
a PA of 75° from the north (the cut is marked in Figure [8j top- 
left panel). The contour levels are from 10 to 90% from the peak 
emission (4.2 Jy beam -1 ) with a step of 10%. The vi sr at 7.7 km 
s _1 and the central position (marked by the cross in Fig.[8]top- 
left panel) is marked by vertical and horizontal line. 



Applying the method of ICabrit & Bertoutl (Il990l Il992l) we 
derived properties of the outflows such as outflow mass, out- 
flow energy and dynamical age. Since the SMA data suffer 
a lot from the missing flux problem, only the 30m data are 
used in the calculations. These calculations assume that the 
13 CO(2 - l)/ 12 CO(2 - 1) line wing ratio is 0. l dChoi et al.lll993b 
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Fig. 8. S255IR 12 CO(2 - 1) outflow images observed with IRAM 30 m telescope and the SMA. The blue-shifted outflow is shown 
in dashed contours with a velocity-integration regime of [-40,0] km s _1 , and the red one is shown in full contours with a velocity- 
integration regime of [16,56] km s~' . The top-left panel presents the SMA data only, the bottom-right presents the IRAM 30 m data 
only. The others are all combined SMA+30 m data, which are inverted with different uv-range selections (shown at the top-right of 
each panel in kA) resulting in different beam sizes shown at the bottom left of each panel. All the 12 CO emission contour levels start 
at 5 cr with a step of 10 cr. The dotted contours are the negative features and the stars mark the position of NIRS 3. In the bottom-left 
panel, the back ground scale is SINFONI H2 line emission. The line in the SMA-only panel shows the pv-cut presented in Fig. [9] 
and the cross marks the central position of the line. In panel (0,120), the SMA 1.3 mm continuum map is over plotted (black full 
contours), and the contours start at 5 cr with a step of 5 cr. In panel (0,70), the same CH3OH (8_i,8 - 7qj)E integrated intensity map 
in Figure |4] is over plotted in black full contours . The cr of the outflow data is shown in Table [7] The (0,0) point in each panel is 
R.A. 06hl2m54.019sDec. +17°59'23.10" (J2000.0). 



Table 7. Outflow data. 



Region 


uv-range 


beam size 


0"red 


cr blue 


Data 




[kA] 


["] 


[mjy beanf 1 ] 


[mJy beam -1 ] 




S255IR 


(8, 173) 


1.5x1.5 


15 


26 


SMA 




(0, 120) 


2.1x1.8 


15 


26 


SMA+30m 




(0, 70) 


3.1x2.5 


18 


32 


SMA+30m 




(0, 45) 


4.0x3.4 


22 


38 


SMA+30m 




(0, 35) 


6.3x5.8 


26 


71 


SMA+30m 




(0, 20) 


9.8x7.4 


38 


86 


SMA+30m 




(0, 16) 


11.1x10.7 


40 


110 


SMA+30m 






11.3x11.3 


0.8Kkms-' 


1.1 Kkms 1 


30m 


S255N 


(8, 173) 


1.8x1.4 


19 


12 


SMA 




(0, 40) 


4.2x3.9 


34 


18 


SMA+30m 






11.3x11.3 


0.5 K km s" 1 


1.1 Kkms 1 


30m 


S255S 


(8, 173) 


1.8x1.4 


14 


18 


SMA 




(0, 40) 


4.2x3.9 


20 


26 


SMA+30m 






11.3x11.3 


0.4Kkms-' 


0.5 Kkms- 1 


30m 



lLevreauldfl988l iBeuther et alll2002h . The results are shown in 
Table [8] The dynamical age calculation depends on the tangent 
of the outflow's inclination angle (with 0° being in the plane of 
the sky). Since we cannot get the inclination angle information 
of these sources, we assumed an inclination angle of 45° for all 



the dynamical age calculations. For S255N and S255S, the out- 
flows are not clearly defined in the 30m-only maps (Figures [TOl 
and [12]). Therefore, we used with the SMA+30m map to deter- 
mine the outflow sizes (see right panels of Figures [TOl and [T2l. 
For S255S, since we just detected only one outflow lobe asso- 
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Fig. 10. S255N 12 CO(2 - 1) outflow images observed with IRAM 30 m telescope and the SMA. The blue-shifted outflow is shown 
in dashed contours with a velocity-integration regime of [-44,0] km s _1 , and the red one is shown in full contours with a velocity- 
integration regime of [16,48] km s _1 . The left panel presents the SMA data only, the right presents the IRAM 30 m data only. The 
middle one is combined SMA+30 m data, the uv-range selection is shown at the top-right of each panel in kA. All the CO emission 
contour levels start at 5 cr with a step of 6 cr. The dotted contours are the negative features and the asterisks mark the position of 
co ntinuum peaks. In the middle panel, the crosses mark t he position of the Class I 44 GHz (7o - 6i) A + methanol masers detected 
bv lKurtz et al.l d2004l) and th e triangle is the water m aser dCyganowski et al.ll2007l) . In the top left panel, the back ground grey scale 
is VLA Q-band continuum fcvganows ki et al.l l2~007) and the line shows the pv-cut presented in Fig.QT] which follows the masers 
direction and crosses the continuum peak. The arrow in the 30m-only panel marks the outflow size we used to calculate the outflow 
physical parameters. The cr of the outflow data is shown in Table [7] The (0,0) point in each panel is R.A. 06hl2m53.669s Dec. 
+ 18°00'26.90"(J2000.0). 



ciated with each continuum source, we just extended the only 
one lobe to the opposite side of the source to get the proper size 
and mass of the outflow. The arrows in the 30m-only panels in 
Fi gures [TOl and[T2l mark the outflow size we use. 

3.4. Rotational structures 

S255IR Figure [13] shows the velocity moment maps of 
HCOOCH3, CH 3 CN (k = 2) and C ls O(2 - 1). In the velocity 
maps of HCOOCH3 and CH3CN (k = 2), we see a clear velocity 
gradient perpendicular to the outflow axis, which indicates the 
existence of a rotational structure. The C 18 0(2 - 1) velocity map 
shows a little different picture compared to the other ones. C 18 
traces a more diffuse gas compared to the other two molecules. 
The C 18 (D velocity map also shows a big velocity difference be- 
tween S255IR-SMA1 and S255IR-SMA2. 

Figure [14] shows the position-velocity diagrams of the 
HCOOCH3 and C 18 0(2- 1) emission. The cuts, which are shown 
in Figure[T3] go through the peak of the dust continuum and have 
position angles perpendicular to the direction of the outflow pre- 
sented in Figure[8] The pv-diagram of HCOOCH3 shows that the 
rotational structure is not in Keplerian motion, hence maybe it is 
jus t a rotating and infalling core similar to the toroids described 
by ICesaroni et alj {2007). The approximately 2.3" diameter of 
the rotational structure corresponds to 3 700 AU at the given 
distance of 1.59 kpc. The pv-diagram of CH3CN (k = 2) has a 
similar structure to the one of HCOOCH3, thus we do not show 
it. 

However, the picture is different for C ls O, the pv-diagram 
shows a Keplerian-like rotation structure. The full line in the 
C 18 panel shows a Keplerian rotation curve with a central mass 



of 28 M Q , which is the mass of the whole continuum structure 
in S255IR covering both mm peaks. For a more detailed discus- 
sion, see Sec. 14.31 

S255N Figure [15] shows the velocity moment maps of 
CH 3 OH(8-i, 8 - 7 , 7 ), CH3OH (8 0; 8 - 7i >6 ) and C ls O(2 - 1). The 
C 18 velocity map shows a complicated velocity structure, the 
velocity gradient is neither aligned with nor orthogonal to the 
outflow orientation. Because C 18 is an isotopologue of 12 CO, 
it may be influenced by infall and the outflow. Methanol is well 
known as a molecu le that traces cores, shocks and masers in star 
formation regions (lJ0rgensen et al.l [2004b iBeuther et al.l l2005b: 
ISoboley et al.ll2007[). It has als o been reported as a low mass disk 
tracer ([Goldsmith et al. 1999). In the velocity maps of the two 
methanol transitions (left and middle panel of Figure [TBI ), we 
see a velocity gradient perpendicular to the outflow axis, which 
indicates the existence of a rotational structure. 

Figure[T6lshows the position-velocity diagram of the (80,8 - 
7^6) methanol line emission, and the pv-diagram of the other 
methanol transition shows the similar velocity structure. The 
cuts go through the peak of the dust continuum and have posi- 
tion angles perpendicular to the direction of the outflow (Figure 
ITST ). The pv-diagram shows that the rotational structure is not 
Keplerian. The structure has a size of 4.4" corresponding to 7 
000 AU at the given distance of 1.59 kpc. It has a narrow ve- 
locity range of 3 km s _1 , and may likely be a lar ge rotating and 
infalli ng core similar to the toroids described by Cesaron i et al.l 
J200l . 

There are two low-velocity components at the northeast and 
southwest of the continuum peak in the methanol velocity maps 
which coincide with the outflow. They may be caused by the 
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Fig. 12. S255S 12 CO(2 - 1) outflow images observed with IRAM 30 m telescope and the SMA. The blue-shifted outflow is shown 
in dashed contours with a velocity-integration regime of [-9,1] km s _1 , and the red one is shown in full contours with a velocity- 
integration regime of [14,32] km s _1 . The left panel presents the SMA data only, the right presents the IRAM 30 m data only. 
The middle one is combined SMA+30 m data, and the uv-range selection is shown at the top-right of each panel in kA. All the 
12 CO(2 - 1) emission contour levels start at 5 cr with a step of 4 <x, except the one in IRAM 30 m only image which starts at 5 cr 
with a step of 2 cr. The dotted contours are the negative features and the stars mark the positions of continuum peak S255S-SMA1 
and S255S-SMA2. The arrows in the 30m-only panel mark the outflow sizes we used to calculate the outflow physical parameters. 
The cr of the outflow data is shown in Table [7] The (0,0) point is each panel corresponds to position R.A. 06hl2m56.58s Dec. 
+17°58'32.80" (J2000.0). 

Table 8. Outflow parameters. 



Sources 


Af, 


P 


E 


Size 


t 


M out 


F m 




[M ] 


[Melons- 1 ] 


[erg] 


[pc] 


[yr] 


[M /yr] 


[M /km/s/yr] 


S255IR-SMA1 


2.9 


137 


6.5xl0 46 


0.7 


7.4xl0 3 


4xl0" 4 


2xl0" 2 


S255N-SMA1 


1.0 


43 


1.9xl0 46 


0.2 


2.4xl0 3 


4x10"* 


2x1 0~ 2 


S255S-SMA1 


0.02 


0.4 


8xl0 43 


0.1 


2.5xl0 3 


8xl0~ 6 


2xl0~ 4 


S255S-SMA2 


0.03 


0.8 


2X10 44 


0.07 


1.7xl0 3 


2xl0~ 5 


5xl0- 4 



Entries include total outflow mass M t , momentum p, energy E, size, outflow dynamical age t, outflow rate M out , mechanical force F, 



shock heating emission, and the line prof ile is consist with 
both thermal emiss i on and maser emission (|Sutton et alll 2004 ; 
ISobolev et aT]|2005b iKalenskii et al.ll2"002tlSlvsh et al.ll2002l) . 



3.5. Temperature from CH 3 CN(12 k - ll k ) in S255IR 

Since we detected 7 lines of the CH3CN(12^ - 1 U) fc-ladder with 
k = 0...6 in S255IR, we can utilize the varying excitation lev- 
els of the lines with lower level energies Ei ower /k between 58 
to 315 K (Table [3} to estimate a temperature for the central gas 
core. Figure [TTI shows the observed CH3CN(12* - 11*) spec- 
trum toward the continuum peak S255IR-SMA1 with only the 
compact configuration data. We did a simple Gaussian fitting of 
the spectrum and plot the level populations Nj^ we calculated 
from the fitting resu lt in Figure [18] wit h the assumption of opti- 
cally thin emission (Zhang et al. 1998). The linear fitting result 
of the lower five levels is also plotted in Figure[l8] It is clear that 
we can not fit the whole spectrum with one single temperature, 
which reveals a temperature gradient of the source and not opti- 
cal thin emission. We modeled this spectrum in the local thermo- 
dynamic equilibrium using the XCLASS software developed by 
Peter Schilke (private communication). This software package 



uses the line catalogs fro m JPL and CDMS dPovnter & Picked 
119851: iMiiller et al.ll200ll) . The model spectrum with a temper- 
ature of 150 K is shown in Figure [17] in dotted line. The main 
difference between the model spectrum and the observed spec- 
trum is that the model one is optically thin whereas the lower line 
intensity of the observed k — 3 line indicates moderate optical 
depth of the CH 3 CN lines. 

3.6. SINFONI results 

3.6.1 . SINFONI line and maps 

The SINFONI observations are centered on the mm peak 
S255IR-SMA1 covering most part of the near-infrared cluster. 
Figure [19] shows the 3 color composite of 3 line maps (Red: 
Bry line emission, green: H2 (2.12 yum) line emission and blue: 
Fe II (1.64 fim) line emission). The two massive young stars 
(NIRS 1 and NIRS 3) are in the center of the field, surrounded 
with the cluster of about 120 stars down to K =17 mag. For 
all the point sources marked in Figure [24] a SINFONI H- and 
/T-band spectrum is available and a spectral classification is ob- 
tained of the brighter members (Sec. 13.6.2b . The PDR and jet 
like emission is traced by the molecular H2 emission (2.12 //m, 
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Fig. 13. HCOOCH3, CH 3 CN (k = 2) and C ls O(2 - 1) velocity (1st) moment maps overlaid with the SMA 1.3 mm dust continuum 
of S255IR. The contours start at 5o~ and increase with a step of 10cr in all panels (<x=1.7 mJy beam 1 ). The lines in each panel show 
the pv-diagram cut presented in Fig. [14] All moment maps were clipped at the five sigma level of the respective line channel map. 
The synthesized beam is shown in the lower left corner of each plot. The (0,0) point in each panel is R.A. 06hl2m54.019s Dec. 
+ 17°59'23.10"(J2000.0). 

S255IR 




Velocity [km/s] 

Fig. 14. Position-velocity diagrams derived for the cuts along the observed velocity gradient in Fig. [13] The offset refers to the 
distance along the cut from the dust continuum peak. The contour levels are from 3cr with a step of lcr in both panels (lcr=60 mJy 
beam -1 ). The full line in the C ls O panel shows a Keplerian rotation curve with a central mass of 28 M Q , which is the mass of the 
whole continuum structure in S255IR. The negative features are shown in dashed lines. 



green). And to the north of the S255IR-SMA1 (the asterisk in 
the Figure[T9), two point sources (sources #17 and #18 in Figure 
l24l) show stron g Bry, which may in dicate the existence of accret- 
ing signature dMuzerolle et alj!998h . The two strong Fe II (blue) 
features around S255IR-SMA 1 show destructive shock (J-type, 
iHollenbac h & McKee ( 1989)) emission and also follow the di- 
rection of the outflow and the H2 jets, which indicate S255IR- 
SMA1 is the energetic driving source of the jets. 

To identify the excitation mechanism of the H2 emission 
(the arc which is marked in Figure [19] by the red contour and 
the jet like emission regions (a) and (b)), we extracted the 
spectra of these three regions and construct the excitation di- 
agrams for these regions. Figure [20] shows the excitation dia- 
grams of them. In these diagrams, the measured column den- 



sities of lines are plotted against th e energy of the upper level 
(see lMartfn-Hernandez et al] d2008l) for the detailed description 
of this diagram). Th e total column dens ities were calculated us- 
ing the description of Zhang et al.l d!998h . Different symbols rep- 
resent different vibrational levels (Fig. [20b. The arc region (Fig. 
[T91 has the most lines detected as it covers the largest area. For 
the jet-like emission region (a) and (b), the weaker lines seen 
in the spectrum of the arc region are not detected and 3cr upper 
limits are given instead. The solid line in the diagrams is a single 
temperature fit to all the data points, while the dashed line in the 
diagram of the arc region presents a linear fit to only the 1-0 
S lines. For the arc region it is clear that the column densities 
are not represented by a single temperature gas, the line fluxes 
of the 2-1 and 3-2 vibrational levels are higher than expected 
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Fig. 15. CH 3 OH(8_i ;8 - 7 , 7 ), CH 3 OH (8 ,g - 7i j6 ) and C I8 0(2 - 1) velocity (1st) moment maps overlaid with the SMA 1.3 mm 
dust continuum of S255N. The contours start at 5<x and increase with a step of 10cr (lcr=1.6 mJy beam -1 ). All moment maps were 
clipped at the five sigma level of the respective line's channel map. The lines in each panel show the pv-diagram cut presented in 
Fig. [14] The synthesized beam is shown in the lower left corner of each plot. The (0,0) point in each panel is RA. 06hl2m53.669s 
Dec. +18°00'26.90" (J2000.0). 



Table 9. Physical properties of the H2 gas. 



Region 


1 rot 




N(H 2 ) 




IK] 




[cm -2 ] 


the arc 


2416 ± 64 


2.6 


±0.2 x 10 17 


a 


1655 ± 62 


5.8 


± 1.7 x 10 17 


b 


1091 ±91 


1.9 


± 1.2 x 10 18 



from the 1-0 line fluxes. A likely excitation mechanism of the 
H2 gas in the arc is fluorescence by non-ionizing UV photons 
(Davi s et al.ll2003l 

The excitation diagrams of regions (a) and (b) are well- 
represented by a single temperature. However the 3-2 lines are 
not detected in either region, and only one 2-1 line is detected in 
region (a) left most of them only upper limits. We conclude that 
the excitation diagrams of both regions are consistent with shock 
excited emission in an outflow emission. Besides, the bottom left 
panel of Figure[8]shows region (a) and (b) follow the direction of 
the outflow, also their elongated shapes that all suggest outflow 
origin. The linear fits of the excitation diagrams also allow us 
to determine the column density and temperature of the emitting 
gas (Table©. 

3.6.2. SINFONI stellar spectral type classification 

Our SINFONI observations provide an H- and Zf-band spectrum 
of each source with a spectral resolution of R =1500. SO FI /- and 
Zf-band photometry results are taken from iBikl (12004 . Objects 
with Zf <14 have high enough S/N spectra to obtain a reliable 
spectral type. This reduces the sample to 39 sources, however, 
we can only get the spectral type of 16 sources which are listed 
in Table QTJ For many point sources, we could not get a spec- 
tral type, some of them have very flat spectra and basically no 
absorption lines and also no infrared excess, they are likely to 
be gas clumps. Others have very red spectra, which are likely 
dominated by dust emission from the surrounding environment 



or circumstellar disks, therefore, the spectrum of the underlying 
objects might not be visible, such as sources #17, #18 in Table 
fTTI and our infrared sources NIRS1 and NIRS3 which will be 
discussed later. 

The SINFONI //+ZT-band spectra of the three brightest clus- 
ter members (sources #1, #2 and #3) show Bry absorption in the 
Zf-band and BrlO - 14 absorption in the //-band (Figure |2TT>. 
Weaker lines of He I are also visible in H (1.70yum) and K (2.058 
ami). We apply the Z T-band classification scheme for B stars from 
lHanson et al. l (fl996T) which links a Zf -band spectral type to an op- 
tical spectral t ype. F or th e //-band, we used the classification of 
lHanson et all (11998) and lBlum et ail (1 19971) . Table [TOl shows the 
measured Equivalent Widths (EW) of the relevant lines in the H- 
and Zf -band spectra of the three B/A stars. Star #1 has a strong 
Bry absorption but does not show He I (2.11 /mi), therefore is 
classified as B3 V-B7V. Star #2 and #3 show much stronger Bry 
absorption and also do not have He I (2.11 /mi) are therefore 
classified as B8V-A3V. The EWs of the //-band lines are in 
agreement with the Zf-band spectral type. Source #1 and sourc e 
#2 are also associated with UCHll regions (ISnell & Ballvf l986). 
which suggests that they are high to intermedia mass stars, and 
this consists with our spectral type classification results. The re- 
sults are listed in Table fTTI 

Besides the 3 B/A stars we found 13 stars showing absorp- 
tion lines typical of later spectral type (Figure l22l . The most 
prominent lines we used are the CO first overtone absorption 
bands between 2.29 and 2.45 //m, and absorption lines of Mg I 
and other atomic species in the //-band as well as Ca I and Na 
1 in the ZT-band. To classify th e late-t ype s tars we use the refer - 
ence spectra of Cushin g et alj d2005l) and iRavner et alj J2009). 
The atomic lines, such as Mg I and Na I, are used for the temper- 
ature determination, while the CO line s are used to determine the 
luminosity class. See Bi k et alj d2010l) for a detailed description. 

Compared with the reference spectra, the CO (2.29 /mi) ab- 
sorption of our SINFONI spectra is usually deeper than in dwarf 
reference spectra, but not as deep as in the giant spectra. In a few 
cases where a luminosity class IV reference spectrum was avail- 
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Table 10. Equivalent width and /iT-band spectral types of the early type stars in S255IR. 
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Star 


Brll (1.68 /im) 


He I (1.70 /mi) 


Hel(2.11/mi) 


Bry (2.166 /mi) 


AT- Spectral type 


Optical Spectral type 


1 


<0.4 A 


0.6±0.4 A 


<0.4 A 


5.9±0.4 A 


KB4-B7 


B3V-B7V 


2 


8.3+0.4 A 


<0.4 A 


<0.4 A 


11+0.4 A 


KB8-A3 


B8V-A3V 


3 


8.9±0.4 A 


<0.4 A 


<0.4 A 


8.9+0.4 A 


KB8-A3 


B8V-A3V 




20 
Velocity [krn/s] 



-20 



b 
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Fig. 11. Position-velocity diagrams of S255N for the 12 CO(2 - 
1) SMA-only outflow observations with a velocity resolution of 
1 km s _1 . The top panel is the pv-plot in northeast-southwest 
direction with a PA of 54° from the north (the cut is marked as 
the line a in Figure [TO] left panel), and the bottom panel is the 
pv-plot in northwest-southeast direction with a PA of 7° from the 
north (the cut is marked as the line b in Figure[10] left panel). The 
contour levels are from 10 to 90% from the peak emission (4.15 
Jy beam -1 in the top panel and 2.82 Jy beam -1 in the bottom 
panel) with a step of 10%. The v; jr at 8.6 km s -1 and the central 
position (i.e. the S255N-SMA1 position in the top panel and the 
cross point of line a and line b in Figure [K)]top-left panel in the 
bottom panel) are marked by horizontal and vertical line. 
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Fig. 16. Position-velocity diagrams derived for the cuts along the 
observed velocity gradient in Fig. Q3] The offset refers to the 
distance along the cut from the dust continuum peak S255N- 
SMA1. The contour levels are all from 3cr and with a step of lcr 
(lcr=40 mJy beam -1 ). 




2.2065 10 2.207 10 2.2075 10 

Rest Frequency (MHz) 

Fig. 17. CH3CN(12^.- 11^) spectrum toward the mm continuum 
peak S255IR-SMA1. The dotted line shows a model spectrum 
with r rot =150K and Nch 3 cn = 3.5 x 10 14 cm -2 . The k = 6 line is 
excluded in the model because it is blurred by the HNCO(10i i 9 - 
9i,g) line. 



able, the spectrum provided a better match to the observed spec- 
trum. This suggests that our l ate type stars ar e lo w- and interme- 
diate P MS stars, and indeed, Luhmanl d 19991) and lWinston et all 
(2009) find that PMS spectra have a surface gravity intermedi- 
ate between giant and dwarf spectra. If the stars were giants, 
dust veiling could make the K-band CO lines weaker. However, 
the H-band CO and OH lines are also much weaker, while the 
atomic lines have the expected EWs. Therefore, this seems to be 
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150 200 
E„A (K) 



Fig. 1 8. From the G aussian fitting result, we calculated the Nj^ 
dZhang et al.| [l998). and plotted in the figure above. E„ is the 
upper level energy of each transition. The line shows the linear 
fitting of the first five points. 
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Fig. 19. Three-color image created from SINFONI line maps. 
Red: Bry line emission, green: H2 (2.12 yum) line emission and 
blue: Fe II (1.64 fim) line emission. The cross marks the posi- 
tion of S255IR-SMA1. The contour and the dashed lines mark 
the regions of which the excitation diagrams are constructed to 
identify the nature of the emission CFig. [20b. 



a surface gravity effect, suggesting a PMS nature of our late type 
stars. 

The "double blind" procedure was applied during the classi- 
fication of the late type stars to check the accuracy of the clas- 
sification. In this procedure, Y. Wang and A. Bik did their own 
classification separately and independently without knowing re- 
sults of each other, and their results were compared afterwards 
to get the differences which are the errors of the spectral type 
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Fig. 20. The excitation diagram of the whole arc in the south- 
east of the Fig. [19](the top panel) and two jet-like emissions (a) 
and (b). The solid line in all diagrams is a single temperature fit 
to all the lines detected, while for the dash-dotted line in the top 
panel a only the 1 -0 S lines are included in the fit. The excitation 
diagram of the arc suggests that UV fluorescence is the excita- 
tion mechanism in this area. Regions (a) and (b) are most likely 
outflows as their excitation diagrams suggest thermal excitation. 



in Table QT| Our classification results showed an error of 1 to 2 

subclass in spectral type. 

The relation from iKenvon & Hartmannl (1 19951) was used to 
convert the spectral type into effective temperature. However, 
this relation applies only for main sequence stars, for PMS stars 
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a different relation may hold (lHillenbrandlll"997L IWinston et al. 
120091) . ICohen & Kuhil (fl979h show that the temperature of PMS 
stars might be overestimated by values between 500 K (G stars) 
and 200 K (mid-K). We took this source of error into account 
when calculating the errors in the effective temperature. 

With the knowledge of the spectral type, we can estimate the 
extinction from the observed color. As the intrinsic J - K color 
for dwarfs and gi ants can differ as much as 0.4 for late K stars 
dKoornneefjl 1983b . for the PMS stars, we used as the intrinsic 
J - K color the mean of the color for dwarfs and giants. The dif- 
ference between the mean value and the dwarf and giant values 
is used as the error in the intrinsic color. The derived extinctions 
are listed in Table [TT] 

3.6.3. YSOs 

Source #17 and #18 have very red spectra and show strong Bry 
point emission, but the spectra are featureless due to the strong 
dust emission, which indicates that they are very young and with 
ongoing accretion activity (Figure l23l. Source #17 also shows 
strong A"-band CO line emissions (Bik & Thlll2004l) . which in- 
dicates the existence of circumstellar disk. The two infrared 
sources NIRS3 and NIRS1 which coincide with our mm sources 
S255IR-SMA1 and S255IR-SMA3, respectively, also have ex- 
tremely red spectra, only several H2 lines can be seen (Figure 
l23l . Therefor we cannot get the spectral type of these sources. 

3.6.4. Cluster membership and the HR diagram 

For the late type stars, as discussed in Sec. 13.6.21 most of them 
have a surface gravity typical for PMS stars, so they are con- 
sidered to be the cluster members. Source #15 and #16 show a 
luminosity type Giant and Super Giant, they are unlikely to be 
associated with this young cluster and since they show quite high 
extinction, they are considered to be background stars. 

To constrain the mass of the cluster members and the age 
of the whole cluster, we construct a HRD to compare the ob- 
served parameters with PMS evolutionary tracks. The derived 
extinctions allow us to de-redden the A"-band magnitude, con- 
vert to absolute magnitude with the derived temperature, and 
plot the points in a A" vs. log(T e ff) diagram (Fig. [25] top). We 
exclude the background stars #15 and #16 to enable a compar- 
ison with the isochrones. The over-plotted dashed line is the 2 
Myr m ain sequence isochrone taken from iLejeune & Schaered 
(2001). The solid line s in the left panel rep resent the evolution- 
ary tracks t aken f rom iDa Rio et"al] (2009) using the models of 
ISiess et all (120001) for 8 different masses: 5, 4, 3, 2, 1.5, 1, 0.5, 
0.4 M . Comparison of the location of the PMS stars with these 
evolutionary tracks yields an approximate mass varying from ~ 
5 M Q for the brightest stars to 0.4 M . 

In the bottom panel of Figure[25] the over-plotted solid lines 
are the PMS isochrones between 0.1 and 10 Myr. The location 
of the stars is consistent with a range of age. Besides source #1, 
#14 and #4 which lie between the 0.5 and 1 Myr isochrones, 
most of the objects span the 1-3 Myr isochrones, with the more 
massive objects closer to the 1 Myr isochrone. For the less mas- 
sive objects, the spread in age is larger, most likely due to higher 
uncertainties in the spectral type determination. The location of 
the PMS stars suggests an age of 2 ± 1 Myr. 

Comparison of the location of the PMS stars in the HRD 
with those of the Herbig AeBe stars (Ivan den Ancker et all 1 998) 
shows that the late spectral type PMS stars are younger than the 
Herbig AeBe stars and will evolve from their present G- and K 



spectral type to late B, A or early F spectral type when they be- 
come main sequence stars. The three relative early spectral type 
sources (#1, #2 and #3) already evolved to late B early A spec- 
tral type. They are in the transit phase between our late spectral 
type sources and the main sequence. 








R.A. [J2000) 

Fig. 24. The stars listed in Table [TT] marked on Bry 
line+continuum map overlaid with H2 line emission contours. 



4. Discussion 

4.1. Different evolutionary stages and triggered star 
formation? 

The SMA and IRAM 30m data together reveal three massive 
star formation regions with different evolutionary stages. The 
SCUBA 850 /im image (Figure |TJ presents three continuum 
peaks in the whole S255 complex region, one in each sub-region, 
which is S255N, S255IR, S255S, from north to south respec- 
tively. With our high resolution SMA observations, -2500 AU 
at the given distance of 1.59 kpc, we found that all SCUBA 850 
Ami sources fragment in to several cores. 

Mini er et al. 1 120071) suggests S255S to be at a very young 
stage without active star formation, however, our observations 
show outflows associated with the mm sources, Furthermore, the 
virial mass is smaller than the mass obtained from the SCUBA 
850 /im measurement, which implies that the S255S region will 
likely collapse. Since the peak column density is also on the or- 
der of the proposed threshold for hig h-mass star formation of 
1 g cirT 2 dKrumholz & McKedl2^08T) . S255S is at a very early 
stage of ongoing star formation and may form massive stars. 
The single dish continuum properties of S255IR and S255N 
do not show much difference, however, our interferometry and 
NIR observations show us the different properties of S255IR 
and S255N. While the large scale NIR emission shows the ex- 
istence of the NIR cluster in S255IR and not so many NIR 
point sources in S255N, which indicate S255IR is most likely 
in a more evolved stage compared to S255N. For the individ- 
ual sub-cores, more lines are detected at S255IR-SMA1 than 
S255N-SMA1, which indicates a higher temperature and more 
evolved nature of S255IR-SMA1 compared to S255N-SMA1. 
Regarding the kinematic properties, S255N-SMA1 has similar 
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Table 11. Photometric and spectroscopic properties of the detected stars brighter than K — 14 mag. 



Star 


R.A. (J2000) 


Dec. (J2000) 


K 


J - K 


T ff 


Sp. Type 


Lum. class 


A v 




(h in s) 




(mas) 


(mas) 


(K) 






(mas) 


1 


06:12:54.91 


17:59:21.05 


11.26+0.03 


1.7+0.06 


15620+2850 


KB4-B7(B3-B7) 


V 


10.0+0.3 


2 


06:12:55.06 


17:59:28.93 


10.60±0.02 


0.8+0.03 


9800+1590 


KB8-A3(B8-A3) 


V 


4.8+0.2 


3 


06:12:54.68 


17:59:32.82 


11.68±0.04 


2.2+0.09 


9800+1590 


KB8-A3(B8-A3) 


V 


12.8+0.4 


4 


06:12:54.82 


17:59:12.98 


13.39+0.08 


2.3+0.21 


3650+271 


Ml.5+1 


V 


8.2+1.0 


5 


06:12:53.60 


17:59:28.18 


13.40+0.08 


2.9+0.26 


4590+460 


K4+1 


V 


12.3+1.2 


6 


06:12:55.31 


17:59:15.90 


13.17+0.07 


1.7+0.20 


4730+446 


K3+1 


PMS 
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Fig. 21. Normalized SINFONI //-band (left panel) and K-b<md (right panel) spectra of the three B/A stars detected in S255IR. Their 
spectra are characterized by absorption lines of the hydrogen Bracket series and He I (dotted vertical lines). The stars are numbered 
according to Table [TT] 



outflow velocities but a much smaller size of the outflow than 
S255IR-SMA1, which may also suggest a younger evolution- 
ary stage of S255N-SMA1. Among other SMA mm continuum 
cores, S255IR-SMA2, S255N-SMA2 and S255N-SMA3, which 
do not have many line emissions are considered to be at earlier 
evolutionary stages. 

Figure[TJpresents the whole star formation region. The young 
star formation region S255 complex lies between the evolved 



H II regions, S255 and S257. From the morphology of the dust 
structure and the H II regions, it appears that the two H II regions 
pushed gas between them and formed the S255IR dust struc- 
ture and t riggered the star formation, which has als o been sug- 
gested by Bieging etail (l2009t) : iMinier et alJ (120071) . To inquire 
that, we studied the velocity map of our 30m data, however, due 
to the high noise level at the edge of the map, we could not find 
a significant velocity difference between the west and east edge 
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Fig. 22. Normalized SINFONI //-band (left panel) and /if -band (right panel) spectra of the 11 objects showing late-type stellar 
spectral type. And most of them are PMS stars. The vertical lines show the location of some of the absorption lines used for the 
classification of the stars. The stars are numbered according to Table [TTI 



of the dust structure to prove the interaction betwee n the H II 
regions and the dust structure. IChavarria et al. (2008) estimated 
the dynamical age of the two H II regions S255 and S257, which 
is ~ 1.5 x 10 6 yr, similar to our NIR cluster in S255IR. So we 
can not prove the triggering star formation. However, we wit- 
ness outflow interaction between S255IR and S255N (Sec. I3.31 >. 
Since S255N is younger compared to S255IR, S255N may be 
affected by S255IR. 

Further more, the age of the cluster around S255IR we 
obtained from the S INFONI data is 2+1 Myr (Sec. l3A4t . 
IChavarria et al .1 (120081) obtained an age of the larger-scale cluster 



of 1 Myr which is consistent with our result. However, the dy- 
namical age of the high-mass protostar in this region, S255IR- 
SMA1, is ~ 10 4 yr (Table|8]l. Without the knowledge of the out- 
flow inclination angle, we could syst ematically underestimat e 
the age value by a factor of 2 to 5 dCabrit & Bertoutl [l992h . 
Therefore, S255IR-SMA1 should not be dynamically older than 
10 5 yr and S255IR-SMA2 is even younger than S255IR-SMA1. 
The age difference between the NIR cluster and the massive pro- 
tostellar objects indicates that the most massive sources in the 
cluster form last. 
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Fig. 23. SINFONI //-band (left panel) and Zf-band (right panel) spectra of the 4 YSOs. The stars are numbered according to Table 



Mini er et al.l d2007l) suggests a collect-and-collapse and trig- 
gered star formation scenario for S255IR, which is that the B- 
stars (from our observations they are late B stars, i.e. star #1, 
#2, #3 in Figure 124b in S255IR formed through a collect-and- 
collpase process, and triggered the formation of NIRS 1 and 
NIRS 3. The NIR and mm data presented in this paper show 
that the cluster detected in the infrared is likely older than the 
mm sources detected with our SMA data. While the SMA mm 
sources lie in the center of the NIR cluster, the NIR cluster mem- 
bers are more distributed around the mm sources. Based on this 
information, we propose an outside-in star formation scenario, 
which is that the central gas filament collapses under the pres- 
sure of the two H II regions. The collapse of the clump may start 
outside-in under the pressure of the two H II regions, the low 
mass cores need lower density to form, and they formed first at 
the outside region of the clump and collapse to stars. And then 
either the low- to intermediate- mass stars may enhance the in- 
stability of the central high-mass cores and potentially trigger the 
high-mass star formation in S255IR, or the massive cores could 
build up slightly slower and start to collapse afterwards forming 
the most massive stars in this region. Our outflow observations 
show a hint that the energetic outflow from the YSOs in S255IR 
may again have affected the star formation in S255N, but this 
needs further observations to confirm. 

However, the different evolutionary stages of the various re- 
gions and even within each region are quite clear. This suggests 
a sequential star formation. 



4.2. Multiple outflows 

We detected outflows in all three regions (Sec. 13. 3t . In the top 
panels of Figure [8] the outflow emission which is nearby the 
SMA continuum sources does not really follow the NE-SW di- 



rection. We suggest this complicated outflow environment in the 
nearby region of the continuum sources is due to the interaction 
between outflows from NIRS 1 and NIRS 3. The red-shifted lobe 
to the south of the continuum source (bottom panels in Figure[8} 
is most likely driven by NIRS 1, because the north-south bipolar 
reflection nebular which is associated wi t h NIRS 1 follows thi s 
direction to the south dJiang et al.l l2008t ISimpson et al.l 12009). 
and the outflow driven by NIRS 3 should be blue-shifted at this 
direction. 

In Figure [10] the red-shifted gas at the bottom part of line 
(b) seems to be associated with S255N-SMA3, however, another 
blue-shifted feature at the southern part of the map shows up as 
the resolution changes. Figure [7] shows that these components 
are mixed together with the outflows in S255IR, which may in- 
dicate the interaction between these two regions. 

In the youngest region S255S, the velocity of the outflows 
is much smaller compared to the other two regions. Table [2] 
shows that the column density of this region is still relatively 
low, S255S is in very young evolutionary stage and likely just 
start collapsing. This is the likely reason we did not detect very 
energetic outflows like for the other two regions in this region. 

4.3. Disk candidates in S255IR 

NIRS 1, which coincides with S255IR-SM A 3, is reported 
to hav e a polarization disk dJiang et all 120081) . ISimpson et al.l 
(2009) also reported a disk, however, with a slightly different 
interpretation. Because NIRS 1 is relatively more evolved, we 
detected only one unresolved mm continuum source associated 
with only the CO isotopologue lines (Sec. 13. 2\ . NIRS 3, which 
coincides with S255IR-SMA1, is considered to be a high-mass 
protostar based on several signatures (e.g., UCHll region, maser 
emissions, hot core emissions and energetic outflows, see Sec. 
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ture described in iFallscheeretail d2009). This structure is also 
much larger than the jeans length which is ~6 000 AU for a tem- 
perature ~20 K and a density ~ 1 x 10 6 . Our system also exceeds 
the cr iteria to maintain a stable disk described by iRratter et al.l 
d2010h . Therefore, we suggest this structure to be a rotating 
toroid which fragmented into several sources, and form a multi- 
ple system. Similar structures of fragmenting pseudo-disks have 
also b een modeled by iRrumholz et al.l d2009t) and iPeters etakl 
(l20Toh . 

Following these, we propose a star formation scenario for 
this large rotational structure, which is that the rotational elon- 
gated dust structure formed first, and at the central region the 
massive source NIRS 1 started forming. However, this large 
structure is much more massive compared to NIRS 1 and is un- 
stable, then it fragmented into two sources, S255IR-SMA1 and 
S255IR-SMA2. 

The SINFONI source #17 shows strong Bry and CO emis- 
sion in /iT-band (Figure 1231) . which indicates the existence of a 
circumstellar disk. Source #17 might also drive the jet like emis- 
sions at the northern edge of Figure [19] It is interesting to find 
disk signatures toward sources with that different evolutionary 
stages within the same forming cluster. 

4.4. Cores and clumps 

Since the single dish has a much larger beam, which can cover 
all the continuum sources, and 13 CO may be optically thick, the 
observed velocity traces the outer layer of the whole clump and 
shows the mean line-of-sight velocity. The molecular lines we 
used to get the veloc ities for the SMA obs ervations are all dense 
gas and disk tracers dCesaroni et al. I l200l . so what we obtained 
are the velocities of the high mass cores (Tabled- The difference 
between the core velocities and the clump velocities of several 
km s^ 1 (Tabled are very differe nt from the low-m ass core cases 
(e.g. 0.46 km s' 1 in NGC 1336 dWalsh et al.ll2007l) . 0.17 km s 1 
in Oph A |D i Francesco et al. 1 12004 1. This difference confirms 
that massive star formation regions are more turbulent than low- 
mass ones. It further implies stronger peculiar motions of the 
protostars and cores within the clump/cluster gravitational po- 
tential. 



Fig. 25. Top : The extinction corrected K vs. log(T e ff) HRD. 
The /iT-band magnitude has been corrected for the distance 
modulus. Over-plotted with a dash line is th e 2 My r main se- 
quence isochrone from iLeieune & Schaererl (|2001), and with 
the solid lines the p re-main-sequence evolutionary tracks from 
iDa Rio et al.l d2009l) . Bottom : The same data but over-plotted 
with solid lines showing the isochrones from IDa Rio et al.l 
(2009). 



Q]and[3]). We detected a rotational structure coinciding with this 
source perpendicular to the outflow (Figure[T3ll. The HCOOCH3 
position velocity diagram shows that this rotational structure is 
not Keplerian, so it could be a rotating toroid around NIRS 3. 
However, the C ls O pv-diagram shows a much larger Keplerian- 
like velocity structure perpendicular to the outflow (Figure IT4l). 
If this source were at a much further distance or observed with 
worse resolution, this structure could easily be identified as a 
disk. However, with our resolution we resolved two mm source, 
and we know it is not a disk. The C 18 gas size is ~ 2 x 10 4 AU 
and the continuum source has a size of 1.6x1 4 AU at the given 
distance of 1.59xl0 3 pc, which is similar to the rotational struc- 



5. Summary 

Combining multi-wavelength data from mm to NIR wavelength, 
we characterize the different (proto) stellar populations within 
the S255 star formation complex. S255S, S255N and S255IR 
show different dynamical and chemical properties, not only at 
mm wavelength but also at infrared wavelengths, which indi- 
cates they are in different evolutionary stages. With the SMA, 
IRAM 30m and VLT-SINFONI observations, we found outflows 
in all three regions, high velocity collimated ones in S255IR, 
high velocity more confined ones in S255N and lower veloc- 
ity confined ones in S255S. The multiple outflows we found 
in S255IR and S255N suggests a potential interaction between 
these two regions. From the outflow maps, we estimated the dy- 
namical age of the outflow driving sources. Although without the 
information of the inclination angle this dynamical age cou ld be 
underestimated by a factor of 2 to 5 (ICabrit & Bertouti l992). our 
mm sources should not be older than 10 5 yr. 

We detected 25 molecular lines in S255IR, including 7 lines 
of the CH 3 CN(12i -11*) k-ladder with k = 0...6, confirming that 
the hot core nature of S255IR-SMA1. Only 10 molecular lines 
are detected in S255N, including 2 CH 3 CN(12 jt - 11 A .) k-ladder 
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with k = 0, 1, which is indicative of a younger age and colder 
temperature. Besides the 3 CO isotopologue line emissions, only 
diffuse SO emission is detected in S255S. This is consistent with 
different chemical ages. 

High-density tracers like CH 3 CN and HCOOCH 3 show rota- 
tional structures around the most prominent high-mass protostar 
candidates in S255IR and S255N. Furthermore in S255IR, the 
C ls O presents an elongated rotational structure with a Keplerian- 
like velocity gradient perpendicular to the outflow. With a size of 
~ 2 x 10 4 AU, this structure can not be a disk but may be a rota- 
tional toroid which fragments into several sources. 

Near infrared H- and /f-band integral field spectroscopy ob- 
servations were done for S255IR. We identified the excitation 
mechanism of the H2 emission. We derived the spectral type of 
16 stars, and 14 of them are considered to be the cluster mem- 
bers. With the knowledge of the spectral ty pe and the SOFI J- 
and K- band photometry results (Bik 2004) of the cluster mem- 
bers, we constructed an HR diagram to estimate the age of the 
cluster, which is 2+1 M yr . This age is consistent with the result 
Chavarrfa et al. I d2008h obtained. The age difference between the 
low-mass cluster and the massive mm cores indicates different 
stellar populations in the cluster. This also leads to a question, 
do the massive stars in this cluster form last? Our data support 
the idea that massive stars form last. We propose the triggered 
outside-in collapse star formation scenario for the bigger picture 
and the fragmentation scenario for S255IR 

The whole picture of the S255 complex suggests triggered 
star formation, however, but we can not give hard proof at 
present. However, the different evolutionary stages between each 
region and different stellar populations in S255IR are consistent 
with sequential star formation. 
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